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Abstract 

The effective Hamiltonian for the R-parity violating couplings induced nonleptonic 
decay B~ — * K~K° is calculated including leading log QCD corrections running from 
rrig, the mass of the up- type squark, to m&. Experimental upper limit of this decay is 
used to give the bound on the R-parity violating couplings A". 
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1 Introduction 



In the supersymmetric generalizations of the standard model, there is a kind of couplings 
which has no standard model analogue. These novel couplings, called the R-parity (or the 
matter parity) violating couplings are not forbidden by some fundamental principles. 

In a variety of supersymmetric models ||, due to additional horizontal symmetries which 
are broken at some scale between the supersymmetry breaking and the Planck scales, the 
R-parity violating couplings are bounded to be small However, these kinds of theoretical 
bounds are quite model dependent. They depend not only on the choices of the horizontal 
symmetries but also on the concrete power counting rules enforced on the superfields [|J. 
Some model dependent cosmological constraints, e.g. X" << 10~ 7 given by the require- 
ment that at the grand unification scale baryogenesis does not get washed out, can be also 
evaded 0. 

Lacking reliable predictions for these R-parity violating couplings, it becomes important 
to determine or bound them from experiments. In the literature some phenomenological 
analyses have been given based on a variety of experiments, from proton decays to the 
precise measurements of the electroweak interactions at LEP 0] (see also J5|). Because these 
measurements are performed at different scales and the QCD couplings between these scales 
are not small, it is necessary to normalize the bounds at a same scale to extract informations 
at high energy scales (e.g. the grand unification or the Planck scale). This is essential to have 
a better understanding on the origins of these couplings. Note that in some newly performed 
renormalization group analyses || , the running of the R-parity violating couplings has been 
carried out between mz and the grand unification scale. 

Recently, it has been pointed out that nonleptonic decays of the B mesons, such as 
B~ — > K~K° whose decay rate in the standard model is quite small, can also provide useful 
bound on the baryon-number nonconserved R-parity violating couplings A"s ||10|| . In the 
present work, we also concentrate on this process. We will calculate the leading log QCD 
corrections to the effective Hamiltonian, sum up the large logs of the form a s \og(mq/mb) 
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(here q is an up-type squark) using the renormalization group equations. We perform the 
hadronic analyses using the perturbative quantum chromodynamics (PQCD) method |TT| , |12 
and then use the current upper limit of the branching ratio to bound A"s. 



2 Effective Hamiltonian with Leading Log Corrections 

The superpotential for the R-parity violating couplings which do not conserve baryon-number 
is 

w = \y.njih. (i) 

where U and D are right-handed superfields and i,j,k are generation indices. These super- 
fields are defined in the basis in which the quarks are in their mass eigenstates. We consider 
only the lightest charge- 1 right-handed squark (denoted as q) and we will return to this point 
later. Following the effective field theory description, q is integrated out at the scale of rriq 
and the results are two 4-quark operators fl2|: 

Ox = (d R ^s R )(s R ^b R ), (2) 
2 = (s R ^s R )(d R Yb R ). (3) 

Then the effective Hamiltonian is written as 

The Wilson coefficients Ci(fi), C*2(/i) are calculated at the scale of q mass as 

C 1 (mg) = l, C 2 (m ? ) = -1. (5) 

The leading log QCD corrections are characterized by the anomalous dimensions for 
these two operators. Calculating the Feynman diagrams in Fig.l, we get the following 
anomalous dimension matrix: 



= 9l 
7 8vr 2 



/ 



1 3 
3 -1 



(6) 



By applying the QCD renormalization group equations 

M^Q(M)=E(T r )«C i (/i), (7) 

we can get the coefficients of the two operators at some low energy scale \x (here fi ~ m&) 
from those at the scale fi = mq. In the case that q is lighter than the top quark, 

CiK) = I^k^ + ^k-))^ 2 ^^^^)-^^))^- 12723 , (8) 

C 2 (m b ) = I^K^ + ^K-))^ 23 -^^^)-^^))^- 12723 , (9) 

where £ = a s {mq) / a s {rrib) ■ Otherwise threshold effects of the top quark should be accounted 
for by two step running: 

with i] = a s (mq)/a s (m t ), and 

Ci(ra 6 ) = ^(C 1 K) + C 2 K))r/ 2 6/23 + ^(C 1 K) 

where ?7 2 = a s (m t )/a s (mb). Note that by taking = 100,300 or 500GeV, the corrected 
Wilson coefficient C\ at the scale of nib is 1.34, 1.45 or 1.50, which shows large QCD correc- 
tions. 



-c 2 K~))?r 12/21 , (io) 

- C 2 (m,))r 12/2 \ (11) 

-C 2 (m t )) % 12/23 , (12) 

-C 2 (m t ))% 12/23 , (13) 



3 PQCD Calculations of the Decay Rate 

Low energy hadronic transitions can be analyzed in the framework of PQCD [|y], [12|]. In 



the present case, the operator 2 which gives the non-factorizable contributions has the 
same Wilson coefficient as that of Oi (up to a minus sign). Thus it is reasonable to use the 
formalism of Simma and Wyler El where no explicit use of the factorization hypothesis is 
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needed. Following this description [T^] we take the interpolating field of the B meson as 

1 /, 



= ^ ^ (j) B (x)^(fi - m B ) 



(14) 



and those of the K and K° mesons as 



(15) 
(16) 



Here I c is an identity in the color space. We have made the approximation = 0. 

The decay amplitudes can be calculated from Fig 2. With the insertions of the operator 
0\, the calculations of the four diagrams give: 



Al = [\dx][dy][dz]Tr 
Jo 

Al = f\dx][dy][dz]Tr 
Jo 

Al = [ [dx][dy][dz]Tr 
Jo 

A 1 , = [\dx)[dy][dz}Tr 
Jo 



(4^) ^ (4^.) ^K-JixPr 



1>B (4 7a ^ S ) ^ K ~ 1 ^ Pr 



Tr 
Tr 



(4™-) j/ Pr 
j s (4 7 ^) 



J l 9 



I 2 ' 



^B ( i—l a 9s \ i>K~ U-^la.9s \ -TlvpR 



Tt^yPfl]-, (17) 

l 9 



where [dx], [dy], and [dz] denote (dxidx 2 ), (dyidy 2 ) and (dzidz 2 ), respectively, and 

lg = yiP- - XjPb, Pb = X 2 p B - lg, Pd = Z 2 Po + lg, Ps = ~Z\Po ~ lg, Ps = VlP- + lg, (18) 

p_ and po being the momenta of the K~ and the K° mesons, respectively. With the insertions 
of the operator 2 , the results are: 



A 2 



Al 



A 2 



[dx] [dy] [dz] Tr 



[dx] [dy] [dz]Tr 



[dx] [dy] [dz]Tr 



i I \ a \ I \ a \ 



j 2 



[\dx][dy][d 
Jo 



z]Tr 



^y7 a ^J ^k-JuPr^kq U—lag s \ ~jp ilPR 

(\a \ i f \ a \ 

( \ a \ ( \ a \ i 



—i 

l 9 

—i 

¥' 

7?' (19) 
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Performing the trace operations in both the spinor and the color spaces, we find that 
the contributions of A\, A\ vanish due to their color structures. The amplitudes A l d and A\ 
are proportional to ran and thus vanish under the approximation of m# = 0. The amplitude 
A\ is the same as A\ up to a color factor. Denoting A a = Ci{m h )A l a + C 2 (m b )^, we get 

A = -2J q d x dydz^CMM^K~{y)Mz) xi{2xi _ xl _ yi){yi _ — y (20) 

8</^ 1 

A b = - dxdydz—^=C 2 {ij)(i)B{x)(t) K -{y)(j) K o{z) — — , (21) 

Jo 3V6 xi{yi — 

A c = dxdydz—^=C 2 (fi)(l)B(x)(j) K -(y)(j) K o(z)— ^ — — , (22) 

3V6 xi{xi - z 1 )(y 1 - x x y 

A d = 0. (23) 
Now the total decay amplitude reads 

A = %%^(A + A b + A c + A d ). (24) 

2m~ 

The analyses given above are independent of the choice of the wave functions. Below 



we will take the wave functions of the B and K mesons as [11, 12 



Mx) = -^j=8(xi-e B ), (25) 

(f> K -(y) = V3fKVi(l-Vi), (26) 

<j) K o(z) = V3f K z 1 (l-z 1 ). (27) 

Now, after integrating over [dx], [dy] and [dz], we get the final amplitudes: 

A a = -2 (4 + 21n±^-ln^ + nr) G, 

A b + A c = f G, (28) 



where 



(29) 



(30) 



The decay width of B — ► K K° is then 

p _ \KbsKds*\ 2 \Ag + Ab + A c \ 2 

64-7rm ~ 

Taking the _B~ meson lifetime as 1.54 x 10~ 12 s, the branching ratio of this decay can be 
easily found. 



4 Result and Conclusion 



The present experimental upper bound for the branching ratio of the decay B~ — > K~K° 
is 5 x 1CT 5 [13| . This gives the constraints for the R-parity violated coupling |Ag 6s Ag ds *| 
depending on the mass of q (see equations (8) - (13) and (28) - (30)). As we have stated 
previously, we have considered only the lightest up-type squark. When the masses of the 
super symmetric partners of the up, the charm and the top quarks are comparable with one 
another, none of their contributions can be neglected. In this case, there exists the possibility 
of cancellations between the contributions from these up-type squarks, and what we get is 

the bound on | Ei=u,c,t C\KbsKds* / ™% I- 

The excluded region in the parameter space is shown in Fig. 3. For the numerical esti- 
mations, we have used parameters as a s (Mz) = 0.117, the hadronic scale fi = lGeV, and the 
decay constants are taken to be fs = 132MeV and /# = 113MeV. The value of e# charac- 
terizes most of the uncertainties in the PQCD calculations. Note that if we take e# = 0.05, 
the bound without QCD corrections coincides with what is given in [[IIJ. Here we also give 
the bound by using e# = 0.065 which is favored by the rare decay of B —>■ K*j 

As a conclusion, the effects of QCD corrections to the effective Hamiltonian induced by 
the R-parity violating couplings are significant. The bound on a particular combination of 
these couplings has been upgraded by including the QCD corrections. 
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Figure Captions: 

Fig. 1 Diagrams in the calculation of the anomalous dimensions of the operators. The 
black blob represents the effective operator 0\ or O2. 

Fig. 1 Diagrams contributing to B~ — > K~K° in PQCD. The black blob represents the 
effective operator 0\ or O2. 

Fig.3 Bounds on |Ag bs A^ s *| from the experimental upper limit of B~ — > K~K° as a 
function of q mass. The upper two lines represent results without QCD corrections, the 
dash-dotted line with e B = 0.05, the dotted line with e B = 0.065. The lower two lines 
represent results with QCD corrections, with the solid one for e B = 0.05 and the dashed one 
for e B = 0.065. 
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